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ABSTRACT 



Context. The Gum 3 1 hubble containing the stellar cluster NGC 3324 is a poorly- studied young region close to the Carina Nebula. 
Aims. We are aiming to characterise the young stellar and protostellar population in and around Gum 31 and to investigate the star- 
formation process in this region. 

Methods. We identify candidate young stellar objects from Spitzer, WISE, and Herschel data. Combining these, we analyse the 
spectral energy distributions of the candidate young stellar objects. With density and temperature maps obtained from Herschel data 
and comparisons to a 'collect and collapse' scenario for the region we are able to further constrain the characteristics of the region as 
a whole. 

Results. 661 candidate young stellar objects are found from WISE data, 91 protostar candidates are detected through Herschel 
observations in a 1.0° x 1.1° area. Most of these objects are found in small clusters or are well aligned with the HII bubble. We also 
identify the sources of Herbig-Haro jets. The infrared morphology of the region suggests that it is part of the larger Carina Nebula 
complex. 

Conclusions. The location of the candidate young stellar objects in the rim of the H II bubble is suggestive of their being triggered 
by a 'collect and collapse' scenario, which agrees well with the observed parameters of the region. Some candidate young stellar 
objects are found in the heads of pillars, which points towards radiative triggering of star formation. Thus, we find evidence that in 
the region different mechanisms of triggered star formation are at work. CoiTecting the number of candidate young stellar objects for 
contamination we find ~ 600 young stellar objects in Gum 31 above our completeness limit of about 1 Mq. Extrapolating the intital 
mass function down to 0.1 Mq, we estimate a total population of ~ 5000 young stars for the region. 

Key words. Stars: formation - Stars: protostars - ISM: jets and outflows - Herbig-Haro objects - ISM: clouds - ISM: bubbles 
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1. Introduction 

The bubble-shaped H II region Gum 3 1 around the young stel- 
lar cluster NGC 3324 is l ocated ^ 1° north-west of the Carina 
Nebula (NGC 3372; see ISmith & Brooksl l2008l for a recent 
review). While numerous observations of the Carina Nebula 
have been performed in the last few years and provided com- 
prehensive information about the stellar populations as well 
as the cloud properties (|Ypnekura et al. 20051 Smith & Brooks 



2007 



2011 



Kramer et al. 2008; Smith et al. 2010a b; Townsl ey et al . 
Preibisch et al. 2011ab,c, 2012; Salatino et al. 20121) the 



Gum 3 1 region has not received much attention. Despite its inter- 
esting morphology and the publicity of HST images of its west- 
em rim (Hubble News Release STScI-2008-34), the HII region 
and its stellar population remain rather poorly studied until to- 
day. This seems to be related to its celestial position: the close- 
ness to the extremely eye-catching Carina Nebula has always 
overshadowed NGC 3324. 

The physical relation between NGC 3324 and the Carina 
Nebula complex (CNC) is still unclear. Recent distance deter- 



* This work is based in part on data collected by Herschel, an ESA 
space observatory with science instruments provided by European-led 
Principal Investigator consortia and with important participation from 
NASA, and on data observed by VISTA (ESO project number 088.C- 
0117(A)), an ESO survey telescope developed by a consortium of 18 
universities in the United Kingdom, led by Queen Maiy, University of 
London. 



minations of Gum 31 and NGC 3324 yielded values of 2.3 kpc 
for the cluster NGC 33 24 (Catalogue of Open Cluster Data; 
iKharchenko et alJ l2005h and 2.5 ± 0.3 kpc for the clouds in 
and around Gum31 ( Barnes et al.ll2010l) . This implies that the 
NGC 3324/Gum31 region is l ocated at the same distance as the 
Carina Nebula (« 2.3 kpc; see lSmithll2"006d) . 

The diameter of the Gum 3 1 Hllregionis ~ 15' (lOpc). This 
H II region in turn is surrounde d by an expanding H I shell which 
encloses it almost completely dCappa et al.. . 20081). 

The available information about the stellar population of 
NGC 3324 is restricted to the three brightest stars. The brightest 
star of the cluster is the multiple star HD 92206, which is con- 
stituted by t wo 06.5V stars HD 92206 A and B and 08.5V star 
HD 92206 C dMafz-Apellaniz et alj[2004l) . A further very lumi- 
nous star in the region is the AO supergiant HD 92207 (= V370 
Car, V - 5.49). This object is a very luminous (log(L/ Lq) = 
5.56) massive (Mjnit ai ~ 30 Mq) evolved star (age » 7 + 1 Myr; 



iPrzvbilla et al.ll2006 ) that drives a very str ong stellar wind (M - 
1.3 ■ 10'^ Mr. vr '; iKudritzki et"anil999h. T he m embership of 



HD 92 207 is not entirelv clear: IClarial (1 1977b and ICarraro et al 



('2001') assume it not to be part of the cluster However, Forte 



(i1976l from HD 92207 b eing wrapped in a nebul ar shell asso- 
ciated with Gum 31) and iBaumgardt etal] ( l2000l from proper 
motions) determine it to be a member 

ICarraro et al.l (12001 1) identify 25 candidate members of 
NGC 3324 with optical photometry and suggest that the clus- 
ter is very young (< 2 - 3 Myr). With three O-type stars (M > 



1 



H. Ohlendorf et al.: Discovering young stars in the Gum 31 region with infrared observations 



18 Mq), the field s tar initial mass function (IMF) representation 
bv iKroupal (l2002h suggests that there should be w 1500 low- 
mass (0.1 Mq < M < 2Mq) stars present. This implies that the 
vast majority of the stellar population of NGC 3324 is still un- 
known. 

In addition to the optically visible stellar cluster in the 
HIl region, there is a population of young stars embedded in 
the molecular cloud surrounding the Hll region. The stars in 
this p opulation are only seen in infrared images ([Cappa et all 
2008h or traced by their protostellar jets (ISmith et al.ll2010ah . 



Cappa et al] (l2008h sampled point sources in a region of 20' ra- 
dius centred on NGC 3324 from the IRAS, MSX and 2MASS 
point-source catalogues and identify 12 (IRAS), 9 (MSX) and 26 
(2MASS) candidate young stellar objects (cYSOs) using colour- 
colour criteria. Due to the limited sensitivity and angular resolu- 
tion of these data, the currently known few dozen of embedded 
infrared sources represent only the tip of the iceberg; many more 
embedded young stellar objects (YSOs) must be present in this 
area and waiting to be discovered. 

In this paper, we are aiming to characterise the protostellar 
and young stellar population of the stellar cluster, the surround- 
ing Hll region and its environs. Our study is based on Spitzer, 
WISE, Herschel, and VISTA data, which provide much better 
sensitivity and angular resolution than the previously existing 
data sets. 

We will describe the data sets we used in Sect. |2] In Sect. |3] 
we will describe the general morphology of the infrared clouds 
in and around the Gum 31 region and its implications for our 
work. Detections of cYSOs from the IR data, including the iden- 
tification method, and their spatial distribution will be described 
in Sect. |4] Furthermore, we have analysed point-like sources as 
detected with both Herschel and Spitzer and derived their spec- 
tral energy distributions (SEDs; Sect. |5). We have also identi- 
fied likely sources to two previously detected Herbig-Haro jets 
(Sect.|6]l. Inferences will be discussed in Sect.|7] 

2. Observational data 

2.1. Spitzer images and piiotometry 
2.1.1. IRAC 

We are using data taken during the Spitzer cold missi on phase, in 
July 2 008, with the InfraRed Array Camera (IRAC; iFazio et al.l 
|2004|) ("Galactic Structure and Star Formation in Vela-Carina" 
programme; PI: Steven R. Majewski, Prog-ID: 40791); retrieved 
through the Spitzer Heritage ArchiveQ The Basic Calibrated 
Data were assembled into mosaics using version 18.4.9 of 
the MOsaicker and Point s ource Extractor package (MOPEX; 
iMakoyoz & Marleaul l2005h provided by the Spitzer Science 
Centr43 (SSC). From the complete data set recorded in the ob- 
servation programme, we use only a 1 .0° x 1 . 1 ° area correspond- 
ing to the Gum 31 region and a small regio n to the west of it. 
This c omplements the area as described by iGaczkowski et al.l 
(1201 3h to complete the area surveyed with our Herschel obser- 
vation programme (cf. Sect. 12.2) . 

Using the Astronomical Point source Extractor (APEX) 
module of MOPEX we then performed source detection and 
photometry on these IRAC mosaics. Photometry was first carried 
out individually for each image in the stack and subsequently 
combined internally to provide photometry data for the entire 
mosaic. Sets of Point Response Functions (PRFs) available from 



' An overview of the observations used can be found in Table[T] 
^ |http : //ir sa ■ ipac. caltech ■ edu/ 1 



the SSC, chosen appropriately for the time of the observations, 
were employed. Before the mosaics were constructed, outliers 
were removed using the Box Outlier Detection method within 
MOPEX and backgrounds between the tiles were matched using 
the Overlap module. Subsequently, three Spitzer Astronomical 
Observation Requests (AORs) each were merged into a single 
mosaic. 

PRF-fitting photometry was carried out separately for all 
four IRAC bands (3.6 |xm, 4.5 \xm, 5.8 ixm and 8.0 |J.m). This 
photometric information was then combined into a single cat- 
alogue by a simple nearest-neighbour matching algorithm, tak- 
ing the 4.5 ^.m band as reference. (We chose this band as our 
reference band because it clearly is the most sensitive one and 
less likely to be subject to misidentifications. These are com- 
mon especially in the longest-wavelength bands where random 
fluctuations in nebulosities are detected as point-like sources.) 
Thus, the resulting catalogue is rather conservative and excludes 
any object not detected in the 4.5 \im band. On the other hand, it 
minimises contamination. The total number of point-like sources 
detected in at least one band in the study area is 57 828. 

We noticed a slight misalignment between the Spitzer bands 
and so globally shifted the 4.5 ^m, 5.8 ixm, and 8.0 |xm band 
positions by -0.17", -0.1", and -0.13", respectively, in dec- 
lination with regard to the 3.6 i^m band that we found to be 
well aligned with the Two Micron All Sky Survey (2MASS; 
Skrutskie et al. 2006). In the resulting catalogue, the root mean 
square (RMS) of the spatial deviations between the IRAC and 
the 2MASS position is 0.22". 

Due to the very strong spatial inhomogeneity of the cloud 
emission in our maps, the sensitivity for IRAC cannot be pre- 
cisely quantified by a single value. Instead, we characterise it by 
two typical values, the detection limit and the typical complete- 
ness limit. The former is quantified by the faintest sources in 
our sample. For 3.6 |xm, 4.5 ixm, 5.8 la.m, and 8.0 ixm these have 
fluxes of 100 \xiy, 100 |aJy, 330 |aJy, and 280 |aJy, respectively. 

To estimate the completeness limits, we constructed his- 
tograms of the measured source magnitudes (Fig. |3) and esti- 
mate the point where the rise in the source count is no longer 
well-described by a power-law. Although it is not a formal mea- 
sure of completeness, the turnover in source count curves can 
serve as a proxy to show the typical values of the completeness 
limit across the field. In this way we estimate completeness lim- 
its of a! 1 .5 mJy, Si 0.7 mJy, ^ 1.2 mJy and « 1 .6 mJy for 3.6 |im, 
4.5 \im, 5.8 |J.m, and 8.0 [im, respectively. 

Comparing thes e to numerical models of stellar evolution 
dBaraffe et alj|1998h . we find that for an age of ~ 3 Myr the pho- 
tospheric emission of a 1 Mq star is well above the detection 
limit for WISE (Sect. l24l i in the 3.4 |j.m and 4.6 iJ.m bands while 
for IRAC the sensitivity extends down to ~ 0.5 Mq. For younger 
stars this boundary shifts upwards so that with WISE 1 Myr old 
0.5 Mq stars are still detectable while with IRAC we could reach 
down to ~ 0.25 Mq. If we compare the complete ness limits to the 
RADMC model ( iDuUemond & D ominik"2004>) which provides 
calculations of continuum radiative transfer in axisymmetric cir- 
cumstellar dust distributions around a central illuminating star, 
we find that disk masses as low as 0.013 would still be de- 
tectable for a 1 Mq YSO. 



2.1.2. MIPS 

For the detailed analysis in Sect. |5]the IRAC data were comple- 
mented by phot ometry on a Multi band Imaging Photometer for 
Spitzer (MIPS) dRieke et alj2004l) 24 |im map ("Spitzer Follow- 
up of HST Observations of Star Formation in H II Regions" pro- 
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Fig. 1. The Gum 3 1 region and its connec- 
tion to the central Carina Nebula, seen in a 
Herschel RGB image, with the PACS 70 ^m 
image in blue, SPIRE 250 |i.m in green, and 
SPIRE 500 (im in red. The blue box marks the 
1.1° X 1.0° area used for analysis here. The di- 
agonal white line marks the border to the left of 
which we obtained IRAC photometric data. 



Table 1 : Overview of observations with Spitzer and Herschel used in this work. 



Instrument 


AOR 


Observation mode 


Observation time 


Date 


PI 


Spitzer IRAC 


23708160 


IRAC map 


3301s 


19 July 2008 


S. R. Majewski 


Spitzer IRAC 


23699200 


IRAC map 


3057 s 


19 July 2008 


S. R. Majewski 


Spitzer IRAC 


23704320 


IRAC map 


3070 s 


19 July 2008 


S. R. Majewski 


Spitzer IRAC 


23688192 


IRAC map 


3082 s 


19 July 2008 


S. R. Majewski 


Spitzer IRAC 


23701504 


IRAC map 


3086 s 


19 July 2008 


S. R. Majewski 


Spitzer IRAC 


23695360 


IRAC map 


3086 s 


19 July 2008 


S. R. Majewski 


Spitzer IRAC 


23696896 


IRAC map 


3087 s 


19 July 2008 


S. R. Majewski 


Spitzer IRAC 


23706368 


IRAC map 


3088 s 


19 July 2008 


S. R. Majewski 


Spitzer IRAC 


23684352 


IRAC map 


3312s 


19 July 2008 


S. R. Majewski 


Spitzer MIPS 


15054080 


MIPS phot 


9469 s 


12 June 2006 


J. Hester 


Herschel PACS 


1342211615 


SpirePacsParallel 


11889 s 


26 December 2010 


T. Preibisch 


Herschel PACS 


1342211616 


SpirePacsParallel 


12863 s 


26 December 2010 


T. Preibisch 


Herschel SPIRE 


1342211615 


SpirePacsParallel 


11889 s 


26 December 2010 


T. Preibisch 


Herschel SPIRE 


1342211616 


SpirePacsParallel 


12863 s 


26 December 2010 


T. Preibisch 



gramme; PI: Jeff Hester, Prog-ID: 20726) retrieved through the 
Spitzer Heritage Archive. The MIPS images were searched by 
eye for point-like sources coinciding with the previously iden- 
tified IRAC and Herschel (Sect. I2.2l i sources. As this field is 
only sparsely populated with point-like sources in the 24 
band, they were easily found by eye. MIPS has a PSF of FWHM 
6"in the 24 \im band, slightly degrading towards the edges of 
the 5.4' X 5.4' field of view, sampled onto pixels with a size of 
2.55". According to the MIPS Instrument Handbook the 24 \im 
sensitivity is highly dependent on the observed sky region, in an 
ideal case it is expected to be ^ 35iaJy. 

Source detection and point-spread fun ction (PSF) photom e- 
try were then performed with StarFinder (iDiolaiti et al.ll2000l) PI 
The PSF is derived as a median value from four point sources. 

^ |http : //www ■ bo . astro . it/StarFinder/index . htm| 



It is then fit to each identified point source together with a slant- 
ing plane for the background and taking into account the contri- 
bution of bright stars adjacent to the fitting region. In this way 
the sub-sample of MIPS point-like sources defined by IRAC and 
Herschel detection is analysed from the brightest (4570 mJy) to 
the faintest (5 mJy) objects. 

2.2. Herschel far-infrared maps 

We analysed the Gum 31 region using data from our recent 
Herschel far-infrared survey of the CNC. These observations 
were performed in December 2010 (Open time project, PI: 
Thomas Preibisch, Prog-ID: OTl-tpreibis-1), using the parallel 
fast scan mode a t 60"s~'. With simul taneous five-band imag- 
ing with PACS dPoelitsch etalJl2010() at 70 ^tm and 160 ^m 
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Fig. 2: An overview of the Gum 3 1 region in the optical, mid-IR, and far-IR. The ESO Wide Field Imager (WFI) optical image [ESO 
image release esol207] shows the V band in blue, R band in yellow, O III 501 nm emission in green, and Ha in red. 



and SPIRE jGriffin et al.ll2010l) at 250 ^im, 350 ^im, and 500 [xm 
two orthogonal scan maps were obtained to cover an area of 
2.3° X 2.3°. The angular resolutions of the maps are 5", 12", 
18", 25", and 36" for the 70 [im, 160 |xm, 250 [im, 350 [im, and 
500 i^m band, respectively. At a distance of 2.3 kpc this corre- 
sponds to physical scales from 0.06 to 0.4 pc. 

A full description of these o bservations and the su bse- 
quent data processing is given by iPreibisch et al.l (1201 2h and 
iGaczkowski et al.l (l2013h . Detection and photometry of point- 
like sources in Hersc hel bands were carried out with CUTEX 
(iMolinari et al.ll20Tll also used for the Hi-GAL survey), a soft- 
ware package developed especially for maps with complex back- 
ground. CUTEX calculates the second-order derivatives of the 
signal map in four directions (x, y and their diagonals), allowing 
the identification of point-like sources by their steep brightness 
gradients. A more detailed descriptio n of the photometry process 
is given bv IGaczkowski et aT] (l2013h . 



As any source catalogue based on maps with strong 
and highly spatially inhomogeneous background emission, our 
Herschel source lists will miss some faint sources and at the 
same time contain a small number of spurious detections. To 
improve the reliability of the source catalogue, we restricted 
our analysis to those sources that are detected (independently) 
in at least two Herschel maps and for the SED construction 
in Sect. |5] even to those detected in at least three bands. For 
this, the point-like sources detected in each individual Herschel 
ba nd were match e d. Th e matching process was described by 
iGaczkowski et al.l ( l2013h . In the area, we detect 91 point-like 
sources in at least two bands and 59 in at least three bands. 
For each of them we checked coincidence with IRAC-identified 
point-like sources. We excluded all cases where either the spatial 
coincidence of Herschel and IRAC source was not clear or where 
more than one IRAC source appeared as a possible counterpart. 



4 



H. Ohlendorf et al.: Discovering young stars 




[5^8] 



Fig. 3: Histograms of the measured source 

This results in the 16 sources identified both in Herschel and in 
IRAC wavelengths that are subject to the analysis in Sect. |5] 

The detection limit can be approximated by the fluxes of the 
faintest detected sources; this is in the range of a: IJy to ^ 2 Jy in 
our maps. An estimate for the 'typical' completeness limit (i. e. 
the limit above which we can expect most objects in the survey 
area to be detected) can be derived from the histograms of the 
fluxes similar to those described in Sect. l2.1.TI The correspond- 
ing limits are at ~ lOJy, ~ 15 Jy, ~ 10 Jy, ~ 10 Jy, and ~ 6 Jy for 
70 i^m, 160 i^m, 250 |J.m, 350 |J.m, and 500 |J.m, respectively. 

2.3. VISTA near-IR images 

An //-band image of the area around Gum 3 1 was ob tained with 
the VISTA InfraRed CAMera (VIRCAM) dPalton et al. 2006) 
in the night of 15 January 2012 as the first observation of our 
Visible and Infrar e d Survey Telescope for Astro nomy (VISTA) 
(llrwinet al.ll2004t lEmerson & SuflierlandlbOlOl) survey of the 
Carina Nebula complex (ESO project number 088. C-01 17(A)). 
VISTA is a 4-m class wide field survey telescope that provides a 
1.3° X 1.0° field of view. The near-infrared camera VIRCAM 
consists of an array of sixteen individual 2048 x 2048 pixel 
Raytheon VIRGO IR detectors, providing more than 67 mil- 
lion pixels with a nominal pixel size of 0.339"on the sky and 
sensitive in a wavelength range of 0.85 - 2.4 |^m. Since the six- 
teen chips are non-contiguous they produce a set of sixteen non- 
contiguous images called a 'pawprint' . For a contiguous sky cov- 



in the Gum 31 region with infrared observations 

1 0000.0 E ' ' ' I ' ' ' I ' ' ' T 




8 10 12 'A 

[4.5] 



1000.0 r 




8 10 12 

[8.0] 



magnitudes in our Spitzer IRAC catalogue. 

erage six pawprints, offset in x- and y-direction, are combined. 
The resulting 'tile' covers an area of 1 .5° x 1.2° on the sky. 

For each of the six pawprints, at five jitter positions 27 expo- 
sures with an integration time of 2 s each have been obtained. As 
it turned out that the seeing conditions (with an average FWHM 
of X! 1 .7") did not meet the pre-specified quality criteria, the ob- 
servations were terminated after the completion of this observ- 
ing block. The observations in the other filters and for the rest 
of our Carina Nebula mosaic positions were successfully com- 
pleted a few months later and processed by the VISTA Data Flow 
System at the Cambridge Astronomy Survey Unit. However, at 
the time of writing the photometric calibration of these data is 
still in progress. Therefore, we use only the image data for a 
(preliminary) scientific analysis here, but no photometric values. 

A preliminary photometric calibration showed that objects 
as faint as H ^ 18.5 are clearly detectable in our VISTA image. 
This is about four magnitudes deeper than the nominal 2MASS 
Survey completeness l imit for crowded loca tions near the galac- 
tic plane of H ^ 14.5 dSkrutskie et al.ll2006l) . 

2.4. WISE 

We used catalogue data from the Wid e-field Infrared Survey 
Explore r (WISE; IWright et al.l 1201 Ol) Afl-Sky Data Release 
(Cutri & al1l2012h . These data were taken during the WISE cold 
mission phase from January to August 2010. We used the stan- 
dard aperture-corrected magnitudes obtained with apertures of 
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8.25"and corrections of 0.222 m ag, 0.280 mag an d 0.665 mag 
for 3.4 ^m, 4.6 [im, and 12 ^m (ICutri et al.ll2012h . The WISE 
observations have an angular resolution of 6.1", 6.4", 6.5", 
and 12.0" for 3.4 n .m, 4.6 \xm, 12 \xm, and 22 [a.m, respectively 
dWright et al.ll20Tq) . In our survey area the catalogue contains 
20 739 point sources with detection in at least one band. Source 
confusion in the catalogue should not constitute a problem in the 
analysis here as we required a matching radius of 0.6" to match 
WISE to IRAC detections, the same as we used for the IRAC 
inter-band matching. On this scale, we were able to distinguish 
clearly between nearest and second-nearest neighbours. 

In way parallel to the way desciibed in Sect. 12.1.11 we es- 
timate the detection limits to be 90 \jJy, 70 \jJy and 700 jxly 
for 3.4 [a.m, 4.6 [a.m, and 12 |j.m (for the analysis described in 
Sect. 14.21 we do not use the 22 |i.m band). Analogously, we esti- 
mate completeness limits of 4 mJy, 3 mJy and 24 mJy for 3.4 \im, 
4.6 i^m, and 12 |.im, respectively. 



3. Morphology 

Figure[T]shows the far-IR morphology of the Gum 3 1 bubble and 
its environs; the cavity of the HII region is clearly delineated. 
Figure |2] contrasts optical and infrared images of the Gum 3 1 
bubble. Near the centre of the bubble a cluster of stars is just 
discernible in the IRAC image. This is NGC 3324. In the MIPS 
24 [a.m and PACS 70 jam images the warm dust surrounding the 
stars of NGC 3324 forms an arc-like structure that follows the 
shape of the inner bubble wall. 

A number of pillar-like structures extend from the edge of 
the bubble into its inner part, especially from its southern rim. 
Notably, some but not all optical pillars coincide with those seen 
in the infrared. The photodissociation regions that very sharply 
delineate the edge of the bubble are well-observable in the IRAC 
5.8 |im image where their fluorescence under the influence of 
UV radiation can be seen, but also in the MIPS 24 ^m image 
that shows the emission from the dust grains within. 

Figure |4a] shows a colour temperature map for the larger 
Gum 31 region, extending downwards into the central Carina 
Nebula complex. It was derived from the rati o of the PACS 
70 \xm and 160 \im emission as detailed by iPreibisch et al.l 
(I2OI2I) . It shows the temperatures for the H II region being com- 
paratively high (w 30 - 40 K), while those of the surrounding 
clouds are considerably cooler down to 20 K. The warm dust 
in NGC 3324 is clearly seen as an outstandingly blue (x; 40 K) 
patch inside the bubble. 

The G286.2lH-0.17 cluster (cf. Sect. l4A4] i also stands out. It 
is about 8 K warmer than its immediate surroundings. It is also 
much denser, as can be seen from Fig. |4b]which shows the hy- 
drogen colu mn density derived fro m the colour temperatures as 
described bv lPreibisch et al.l(l2012h . We see that the column den- 
sity inside the bubble is relatively low compared to its surround- 
ings at a few 10^° cm"^ and rises steeply by more than one order 
of magnitude at the ionisation front. Other than G286.21+0.17, 
we see some more dense clumps scattered within the bubble rim 
and beyond it. A notable feature is the cluster G286. 38-0.26 (cf. 
Sect. I4.4.5] l which has a column density of around 3 ■ 10^^ cm"^. 
In sharp contrast to the unusually warm clump G286. 21-1-0. 17, 
in Fig. |4a]G286.38-0.26 is seen to be cooler than its surround- 
ings, down to 20 K. The integrated cloud mass derived by 
IPreibisch et al.l is 186 700 Mq. They derive this value by inte- 
grating over the same column density map that we show here, 
but for a region around Gum 31 slightly differently defined than 
the one here. 



In the longer IRAC-wavelength emission and especially in 
the Herschel image, a 'bridge' of filamentary structure can be 
seen to extend from Gum 31 downwards in the direction of 
the central Carina Nebula, forming a connection. The column 
density map (Fig. l4bl i. too, shows that the clouds suiTounding 
Gum 31 are connecte d to the clou ds in the more central parts 
of the Carina Nebula. 'Forte' (1976) remarks that in the deep op- 
tical plates of Lynga & Hansson ( 1972i) a fila mentary struct ure 
connecting the HII region H-31 (^ Gum 31; lHoffleitllT953l) to 
NGC 3372, the c entral Carina Nebula, c an be seen. 

According to lYonekura et al.l (|2005|) . the radial velocities of 
their C"^0 clumps Nos. 2-6, which surround the Gum31 bub- 
ble, range from Vlsr = -20.0kms"' to Vlsr = -24.1 kms"'. 
The C"^0 clumps in the central and northern part of the Carina 
Nebula (Nos. 8 - 12) have radial velocities in the Vlsr = 
-25.8 ... - 19.9 km s"' range. This good agreement of the ra- 
dial velocities suggests that the clouds around Gum 3 1 and the 
Carina Nebula are conn ected and actual l y part of the same com- 
plex. As was argued bv lPreibisch et al.l (1201 2l) and in the intro- 
duction, we will therefore assume a distance of 2.3 kpc towards 
Gum 3 1 , the same as to the Carina Nebula complex. This number 
also agrees well with recent distance determinations with inde- 
pendent means (see Sect.lTJ. 

4. Young stellar objects in the Gum 31 region 

4.1. Identification of YSO candidates 

The selection of cYSOs in the following sections is based on 
the detection of infrared excesses. It follows, therefore, that only 
those YSOs exhibiting excess infrared emission can be detected 
with these methods (cf. also Sect. 14.31 ). Infrared excess emis- 
sion is indicative of circumstellar material. This means that only 
YSOs of Class to Class II are the subject of this analysis while 
Class III objects must remain undetectedQ 

4.1 .1 . Selection by Spitzer colours 

With IRAC data only and following the crite ria established for 
a survey of star- forming region Pismis 24 ([Fang et al.l 1201 2t 
I Allen etal.ll2004l) we were able to identify 304 infrared ex- 
cess sources from plotting the [3.6] - [4.5] colour against the 
[5.8] - [8.0] colour for those 6739 point-like sources detected in 
all four bands. For ide ntification as a cY SO iFang et all (1201 2h 
demand (Eq.mfoUows lAllen et al.ll200l : 

[3.6] - [4.5] > and [5.8] - [8.0] > 0.4 (1) 

[3.6] - [4.5] > 0.67 - ([5.8] - [8.0]) x 0.67 . (2) 

The resultant colour-colour diagram is shown in Fig. |5] 

If we also take into account the 2MASS magnitudes, we can 
sample a much larger portion of sources by applying the excess 
criteria of 'Winston et al.' (20071 using Spitzer extinction as de- 
rived bv .Flaheitv e t al. 2007). These yield another estimate of 
the YSOs in the region from the J - H \s. H - [4.5] colour- 
colour diagram. In this way, 2577 cYSOs are identified (from a 
total sample of 38459 with detections in the three bands). 



YSOs a re commonly classified according to the slope of their SEDs 
l lLadal 19871) , from Class protostars in the main collapse phase through 
Class I YSOs shrouded in envelopes and accreting infalling circumstel- 
lar matter and Class II YSOs, T Tauri stars with disks, to Class III ob- 
jects with little or no circumstellar matter. 
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Fig. 4: Maps derived from Herschel data. They show the cloud structure in and around the Gum 3 1 nebula and the connection to the 
central parts of the Carina Nebula (just outside the lower left edge of the maps). 
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Fig. 5: Colour-colour diagram using all Spitze r IRAC bands. The 
crit eria adopted for YS O identification follow iFang et al.l(l2012l) 
and lAllenetal1(l2004l) . 

Within this sample, it is certain that there are contaminants 
masking as YSOs. AGB stars unfortunately exhibit very simi- 
lar colours to YSOs and so do background star-forming galax- 
ies. lOUveira et al.l (|2009|) in a spectroscopic survey of objects 
discovered by the Spitzer Legacy Program "From Molecular 
Cores to Planet-Forming Disks" ("c2d") in the Serpens molec- 
ular cloud find a contamination of 25% within their subsam- 
ple of the brightest objects (/^s.oum > 3mJy). They attribute 
this high proportion to the closeness of the Serpens cloud to 
the Galactic plane - something which is also true for Gum 3 1 
(^Gum3i ~ 0.2°). We can also estimate the contamination based 
on the criteria of IWinston et al.l (|2007|) for the selection of con- 
taminants from the [5.8] - [8.0] vs. [4.5] - [5.8] and [4.5] - [8.0] 
vs. [3.6] - [5.8] diagrams. This yields an estimated contamina- 
tion in the IRAC-selected sample of < 10%. 

It is noticeable that the distribution of these cYSOs through- 
out Gum 31 and outside it is almost uniform. This is suspicious. 
We thus conclude that the contamination within this sample due 
to fore- and background sources is high and do not use this sam- 
ple for further analysis. 

4.2. Selection by WISE colours 

We used WISE catal o gue d ata to search for infrared excess 
sources. iKoenig et al.l (l20I2h in their study of massive star- 
forming regions developed a set of criteria to identify likely 
YSOs from WISE four-band photometry. We applied these cri- 
teria accordingly. 

Using the [3.4] - [4.6] and [4.6] - [12] colours, partly in 
combination with the [3.4] and [4.6] magnitudes, probable back- 
ground objects are removed from the sample before the YSO se- 
lection. This includes galaxies (very red in [4.6] - [12]), broad- 
line AGNs (of similar colours as YSOs, but distinctly fainter) 
and resolved PAH emission regions (redder than the majority of 
YSOs). From this cleaned sample the IR excess soiirces a re then 
selected by demanding (Eqs.l3]andl4l lKoenig et alJl20I2[) 

[3.4] - [4.6] - 0-1 > 0.25 and [4.6] - [12] - 0-2 > 1.0 , (3) 

where cr is the quadratically added uncertainty of the respective 
magnitudes. Class I sources are a subsample of this defined by 

[3.4] - [4.6] > 1.0 and [4.6] - [12] > 2.0 (4) 



Fig. 6: Colour-colour diagram using three out of four WISE 
bands. The c riteria adopted for YSO identification follow 
iKoenig et alj ( |20I2|) . Probable YSOs are marked in red for 
Class I sources and in green for Class II. 



(the rest are Class II objects). From this analysis we excluded any 
data point from the catalogue that did not have a signal-to-noise 
ratio of 5 or better Figure ||]shows the [3 .4] - [4.6] vs. [4.6] - [12] 
colour-colour diagram, constructed after probable contaminators 
had been removed. The total sample before this removal were 
10 128 sources. After removal of probable contaminators, 6669 
sources remained. These were plotted in Fig. |6] Resultant Class I 
sources are marked in red. Class II in green. The same colour- 
coding is used in Fig. [8] which shows the spatial distribution of 
the cYSOs. This analysis yields 661 cYSOs of which 207 are 
Class I and 454 are Class II. 

The distribution o f WISE-selected cYS Os follows the cloud 
structure as expected. IKoenig et al.l (1201 2l) estimate the remain- 
ing contamination of the cYSO population selected with their 
criteria by galaxies to be ~ lOdeg which in our 1.2 deg field 
leads to an estimate of ~ 12. However, in a region projected onto 
the Galactic disk as the Carina Nebula we expect a high level of 
contamination mainly from background and foreground stars . To 
estimate the total contamination we performed the same analy- 
sis as described above for two circular areas of 30' radius well 
outside the Carina Nebula (i. e. regions that we expect to be rel- 
atively free of YSOs), one centred around aiooo - 10:51:52.5, 
§2000 = -59:14:30, the other around aaooo = 10:25:10.7, 82000 = 
-58:52:00. They were chosen as having a similar proximity to 
the Galactic plane as Gum31 (b - 0.15 and b = -1.18; Gum31: 
b - -0.17). Furthermore, we carefully selected fields that were 
as free of CO and Ha emission as possible and appeared as 
empty as possible in the IRAS 100 |J.m images. 

In those control fields we find a mean cYSO density of 
97 deg"^; with the cYSOs spread homogeneously. In the Gum 3 1 
region we determine a cYSO density of 550 deg"^. This leads us 
to expect a contamination of around 18% (~ 120 cYSOs) for 
the Gum 3 1 cYSO sample. For the IRAC data we were not able 
to conduct a similar comparison as we do not have photometry 
data for the full study area around Gum 3 1 and none for regions 
outside the central nebula. 

We therefore deem the WISE-selected sample to be more re- 
liable and base our conclusion in the following sections predom- 
inantly on this sample, although we will occasionally describe 
a classification with IRAC. In general, though, those should be 
treated as less reliable. 
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4.3. Identification of protostars from Hersctiel data 

With the methods descri bed in Sect. 12.21 and in more detail by 
iGaczkowski et al.l (1201 3h . we obtained a point-source catalogue 
for the Herschel data. Although we consider only 'point-like' 
Herschel sources in the following, it is important to keep in mind 
that the relatively large PSF corresponds to quite large physical 
scales at the 2.3 kpc distance of Gum 31. In the PACS 70|a.m 
map, all objects with an angular [spatial] extent of up to w 5" 
[0.056 pc] are compact enough to appear 'point-like'. For the 
SPIRE 250 ^.m map, this is true for sources of up to ^ 18" 
[0.20 pc]. This shows immediately that (pre-stellar) cloud cores, 
which have typical radii of % 0.1 pc, cannot be (well) resolved 
in the Herschel maps and may appear as compact 'point-like' 
sources. This implies that YSOs in all their evolutionary stages 
can, in principle, appear as point-like sources in our Herschel 
maps. However, the possibility to detect an object in a spe- 
cific stage depends stron gly on its properties: As described by 
IGaczkowski et al.l (l2013h . many pre-stellar cores (> 1 - 2Mq) 
and embedded protostars (> 1 Mq) will be easily detectable, 
while most of the more evolved pre-main sequence stars with 
disks should remain undetected. In this last case the detection 
limit depends on the disk mass; as shown bv IGaczkowski et all 
a 1 Mq YSO is still detectable if it has a disk mass of > 0.5 Mq, 
this sinks to > 0. 1 Mq for a 3 Mq YSO and > 0.01 Mq for a 6 Mq 
YSO. 

This means that while with IRAC we find mainly Class II 
YSOs and a number of Class I YSOs, with Herschel the em- 
phasis is on Class protostars, with some Class I stars. This 
also imp lies that the overlap between both is relatively small (cf. 
Sect.lSt. TRagan et aP (l2012l) presented radiative transfer models 
of starless cores and protostellar cores and investigated the de- 
tectability of these classes of objects. Their models showed that 
the SEDs of starless (i. e. pre-stellar) cores typically peak around 
200 -300 |J.m and drop very steeply towards shorter wavelengths. 
Their model fluxes at 70 |xm (scaled to the distance of the CNC) 
are several orders of magnitudes below our detection limits. 
Protostellar cores, on the other hand, have much stronger fluxes 
at PACS wavelengths. Guided by these results, we can take a de- 
tection at 70 ixm as an indication for the protostellar nature of 
the source, whereas Herschel sources without 70 |xm detection 
would then be pre-stellar cores. In Fig. Qthese two classes are 
ind icated separately. 

IGaczkowski et alJ (l2013h argue that in a sample of Herschel 
point-like objects in the Carina Nebula it is very unlikely for 
contamination by evolved stars or extragalactic objects to occur. 
The same reasoning applies to Gum 31. The photometric data 
for those sour ces that fall within the Gu m 3 1 area can be found 
in the paper of IGaczkowski et aP (l2013l) . 

4.4. Spatial distribution of the cYSOs 

It is notable that both Class I and Class II sources are found pre- 
dominantly within the interior of the Gum 3 1 bubble or along its 
rim. They tend to occur in clusters and i n thei r distribution are 
often correlated with the lYonekura et al.](l2005b molecular cloud 
clumps (see below). When in the following we refer to a cYSO 
with a number preceded by '[CNA2008]' this source was classi- 
fied as a cYSO bv iCappa et al.,(.2008.) . 

4.4.1. NGC3324 

The cluster NGC 3324 appears prominent in the Spitzer images, 
containing ^ 200 point-like sources. Within the cluster itself, we 



identify only a single cYSO with WISE, J103706.9-583710 at 
the western edge of the cluster (Class II). From the [3.6] - [4.5] 
vs. [5.8] - [8.0] diagram (IRAC) we identify half a dozen cYSOs 
coincident with the cluster NGC 3324 and distributed about 40" 
towards the north-east of the cluster centre. The very small frac- 
tion of stars with detectable infrared excess (w 0.5%) supports 
previous age estimates of > 3 Myr for this cluster 



4.4.2. cYSOs in the rim of tine bubble 

Numerous cYSOs are found lined up along the ionisation front 
of the Gum 31 bubble to the west of NGC 3324. A dozen 
cYSOs are found right behind the edge of the ionisation front. 
Among them are J103653.9-583719, J103653.3-583754, and 
J103652.4-583809, three Class I candidates that are found along 
the ridge traced in IRAC and Herschel images at the very edge of 
the bubble, neighbouring NGC 3324. Two of them, J103653.9- 
583719 and J103652.4-583809, will be discussed in Sect. as 
possible sources of Herbig-Haro jets. 

Behind this 'first row' of Class I candidates, there is a 'sec- 
ond row' of five Class II candidates, all lined up about 19" be- 
hind the ionisation front. Further northwards along the rim there 
are four more Class II candidates, one similarly behind the front 
and three right along it. 

In Fig.[8]it is evident that a major part of the cYSOs is located 
along the bubble edges, similar to what is obse rved in compara- 
ble bu bbles associated with H II regions (e. g. iDewangan et al.l 
l2012h . Their distribution follows its circular shape and their 
number sharply goes down outside the cloud structure that ap- 
pears magenta in the image. This is suggestive of triggered star- 
formation in a 'collect and collapse' scenario as described by 
IWhitworthet"al1 (11994 . 

Following analytical models for this scenario, a stellar wind 
bubble produced by three massive stars with spectral types 06.5, 
06.5 and 09.5 could reach a radius of approximately 9 - 1 1 pc in 
1.5 - 2.0 Myr assuming initial cloud densities in a range of 500 
- 1000 cm"^. The radius at which fragmentation occurred would 
then vary between 7 pc and 1 1 pc. To estimate the mechanical lu- 
minosities emitted by the three most m assive stars in NG C 3324, 
we use d values from ISmithI (l2006al) (and the erratum (ISmithI 
l2006bl) ') for stars with the same spectral types (Luminosity class 
V) observed in the Carina Nebula. The ~ 15' diameter of the 
Gum31 HII region corresponds to lOpc at a distance of 2.3 kpc 
which agrees very well with the values derived from the model. 
The fact that we also find numerous embedded cYSOs in and 
near the rim of the bubble is consistent with star-formation ac- 
cording to the 'collect and collapse' model. 



4.4.3. cYSOs in pillars 

In the southern and eastern part of the bubble edge we find a 
number of small pillars extending into the bubble interior They 
can be seen in the optical and infrared in Fig. |2]and a closeup 
of the IR pillars is shown in Fig. |9] Within four of them we 
find cYSOs in their very tips, reminiscent of what is observed in 
the central Carina Nebula (e. g. the South Pillars region). One of 
them, J103806. 6-584002, coincides with a Herschel point-like 
source (cf. Sect. 15.1b and is therefore most probably a protostar 
This suggests that radiative triggering is at work, very similar 
to the processes seen in the South Pillars (ISmith et al.ll2010bt 
IGaczkowski et al.ll2013l) . 
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Fig. 7. Herschel RGB image (red: 500|i.m, 
green: 250 blue: 70 (J-m) with positions of 
all Herschel point-like sources detected in at 
least two bands overlaid. Red crosses show 
the position of protostellar, yellow boxes those 
of prestellar cores. The large yellow rectangle 
marks the borders of the region analysed. 



4.4.4. G286.21+0.17 

The cloud clump G286.21+0.17 ([DBS2003] 127, BYF73) is lo- 
cated about 12' north of the rim of the Gum 31 bubble. The 
structure of t his clump was rece ntly studied in several molec- 
ular lines by iBarnes et al.l (l2010l) . They determined a diameter 
of ~ 0.9 pc, a luminosity of ~ 2 - 3 ■ 10"^ L©, and estimated 
a clump mass of about 20 000 M© (a value that is about forty 
times large r than the previous m ass estimate based on millime- 
tre data by iFaundez et al.l 120041) . Their derived mass infall rate 
of ~ 3.4 ■ 10"^ Mq yr"' would be the highest mass infall rate yet 
seen, if confirmed. These properties make this cloud a particu- 
larly interesting site of possible massive-star formation. 

In om Herschel maps, this clump appears as a very bright and 
prominent compact feature. In Fig.[TTl we compare its morphol- 
ogy in the Herschel far-IR bands to a near-IR image from our 
VISTA data and a Sp itzer mid-IR image. As already discussed by 
iBarnes et al. l (l20To() . a young stellar cluster, surrounded by dif- 
fuse nebulosity, is located immediately north-west of the clump. 
This cluster appears very prominent in Fig. |8]as well and con- 
tains ~ 45 cYSOs. In the centre of the clump itself, the Spitzer 
images show two bright point sources with an angular separa- 
tion of 7.6" (Fig. fTTT i. They can be identified with the 10 ^.m 
point sources J103832.08-581908.9 and J103832.71-581914.8 
that were detected as counterparts of th e MSX/RMS massiv e 
cYSO G286.2086+00. 1694 described bv lMottram etall (|2007|) . 
for which a bolometric luminosity of 7750 Lq was determined. 

The VISTA //-band image shows very faint diffuse nebulosi- 
ties at the location of these two mid-IR sources (see the close-up 
in Fig. [Toll. This is consistent with the idea that these two objects 
represent deeply embedded YSOs in the protostellar evolution- 
ary phase. With WISE the two sources are not resolved but run 
into one that we classify as a Class I candidate. 

The peak of the Herschel PACS 70 \xm emission is centred 
on the mid-IR source J103832. 0-581908. A two-dimensional 
Gaussian fit to the 70 |im emission yields a FWHM size of 
13.6" X 11.3", which is clearly larger than the FWHM size of 



10" X 10" measured for several isolated point-like sources in the 
same map. With an angular distance of 7.6", the emission of the 
two mid-IR sources cannot be resolved in the Herschel maps, but 
the measured direction of the elongation towards J103832.71- 
581914.8 suggests that both contribute to the observed far-IR 
emission. 

A second cYSO is found in the south-east of the clump, at 
the projection of the Une connecting J103832. 08-581908.9 and 
J103832.71-581914.8. It probably corresponds to the IRAC- 
identified source at aaooo - 10:38:33.6, 62000 = -58:19:22. 
There is a third identified cYSO, situated at the very edge of 
the clump, at aaooo = 10:38:35.0, Saooo = -58:18:44. These 
three are the only cYSOs identified in the cluster with WISE, 
while in the immediately adjoining cluster of stars visible in the 
IRAC images a large number of c YSOs is found. There are ~ 45 
overall, with a slight majority of Class II sources over Class I 
sources. If there is any trend in their spatial distribution. Class I 
are found with slight emphasis to the south-east, while Class II 
sources tend to be located towards the north-west. 

Our colour-temperature map (Fig. |4a), constructed from the 
Herschel 70 |.im and 160 \xm maps, shows that the cloud temper- 
atures range from $ 20 K at the edge of the clump to ~ 25 - 30 K 
in the clump centre, and up to 33 K in the nebulosity surrounding 
the stellar cluster north-west of the clump. 

In our Herschel column-density map (Fig. l4b] i. the level 
A^H - 2 ■ 10^^ cm"^ traces the shape of the clump. This agrees 
ve ry well with the morph ology as seen in the 1 .2 mm map shown 
bv iFaijndez et al.l (l2004l) . The peak value of the column density 
is found to be 1.4 • 10^^ cm"^ and corresponds to a visual extinc- 
tion of Av ~ 70 mag. From our column density map we deter- 
mined the mass of the clump by integrating over a 200" x 200" 
(2.23 pc X 2.23 pc) box around the clump and subtracting the 
local background level. This yields a clump mass of 2 105 Mq. 
Th is value is nea r ly five times as large as the 470 Mg derived 
by iFailndez et al.l (l2004j). but a factor of t en smaller than the 
2OOOOM0 estimated bv lBarnes et al.l (l20IOl) . This seems to sug- 
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Fig. 8: Spitzer IRAC RGB image (red: 8.0 yim, green: 4.5 \xm, blue: 3.6 |am) with positions of WISE-identified cYSOs overlaid. 
Class I sources are marked by red diamonds. Class II in green. The white ellipses represent the clusters discussed in Sect. 14.41 with 
the yellow circle marking the multiple star HD 92206 which contains the O stars in NGC 3324. The orange circles mark the cYSO 
clusters discussed in Sect. 14.4.61 The yellow rectangle marks the borders of the region analysed. 



gest that the mass infall rate estimated bv lBames et al 1 (120101) IS 
also too high, and this clump is not as extreme as suspected. 

Its far-IR fluxes as derived from our Herschel data using 
elliptical apertures for photometry with GAIA|f| are 1661 Jy, 
2261 Jy, 1293 Jy, 653 Jy and 297 Jy for 70 ^m, 160 250 ^m, 
350 i^m and 500 |im, respectively, which results in an integrated 
far-IR luminosity (70 |im to 1 .3 mm) of Lim ~ 9000 Lq. 

This is clearly one of the most luminous clumps in the CNC. 
Its mass is rather high, but probably not as high as previously 
suggested. It may form stars with % 10 M©, but is probably not 
massive enough for the formation of high-mass stars with M* ^ 
20 Mq. 



^ |http : // astro . dur . ac . uk/- pdraper/gaia/gaia . html | 



4.4.5. G286.38-0.26 

The Spitzer IRAC images show a prominent dense cluster of 
several dozen stars at the southern edge of the Gum 31 bub- 
ble, which is surrounded by bright diffuse nebulo us emission 
(Fig. [Til l. This cluster is Usted as [DBS2003] 128 bv lDutra et all 
It is spatiall y coincident with the extended (r - 2pc) 
C'^O clumpQNo. 6 jVonekuraet al.ll2005h . 

The nebulosity around the stellar cluster displays a remark- 
able arc-like shape at the eastern edge. Projecting it into a 
full circle, it would have around 42" diameter in IRAC im- 
ages and 62" in Herschel images. The centre of this circle 

* We note that lYonekura et al.l ( 120051) denoted these structures as 
'cores'; however, according to the definition that cores are very com- 
pact clouds (with typical sizes of ~ 0.1 pc or less), out of which individ- 
ual stellar systems form, these clouds are better characterised as chimps 
(i. e. relatively large dense clouds linked to the formation of small stellar 
clusters). 
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Fig. 9: In a Spitzer IRAC 8.0 |^m image, the yellow arrows mark the positions of four remarkable cYSOs found in the very tips of 
pillars. Positions of WlSE-identified cYSOs are overlaid as in Fig. |8] Class 1 sources are marked by red diamonds. Class II sources 
in green. Blue crosses mark Herschel point-like sources detected in at least two bands. 



would be around aaooo - 10:38:03, 62000 - -58:46:19. The 
star J10380461-5846233, cYSO [CNA2008] 21, is found close 
(~ 11") to this central position of the arc. In the Spitzer bands, 
the cYSO shows strongly increasing brightness with wavelength. 
In the MIPS 24 \iva image, it is the brightest point-source in the 
cluster It was detected as a mid-IR source by MSX and is listed 
as G286.3773-00.2563 in the MSX6C catalogue. With WISE 
data we classify it as a Class II cYSO. The star is not detected 
in any of our Herschel far-IR images. Using USNO-B optical 
catalogue data, 2MASS, IRAC, MIPS, and WISE photometry 
an d Herschel upper limi ts we employed the online SED fitter 
bv lRobitaille et alJ(l2007h to construct an SED and thus estimate 
the (proto-) stellar parameters. The stellar mass is estimated to be 
5.8 Mq for the best-fit model, the luminosity a; 238 L©. Within 
the arc three further WISE Class I cYSOs are seen, J103805.8- 
584542, J103758.4-584648 and J103800.7-584654. 

The optically brightest star in th e cluster is HP 303094, for 
which a spectral type A2 is given bv lNesterov et al.l(ll995h . It is 



located about 17" south of J10380461-5846233 and the centre 
of the arc. According to the Pickles & Depagne (2010i) survey of 
all-sky spectrally matched Tycho-2 stars it may be a foreground 
star (distance: 886 pc). 

Strong far-infrared emission from the region of this cluster 
was detected with IRAS (point source IRAS 10361-5830). Our 
Herschel maps resolved this far-IR emission into ten point-like 
sources in the area of the clump. The only Herschel source with 
near-IR counterpart is J103801.4-584641. The IRAC image at 
this point is dominated by strong nebular emission. There is 
extended emission in a confined region a few arcseconds east 
of the //eric/ieZ-identified point-like source, but since it is not 
well-resolved and the identification with the //eric/zeZ-identified 
source is not unambiguous, we do not include it in the sam- 
ple studied in Sect. |5] The two //erac/ieZ-identified point-like 
sources north-east of the arc-like nebula that are also detected 
as bright sources in the MIPS maps are J103810.2-584527 and 
J103807. 2-58451 1, both have no clear near-IR counterparts. A 
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Fig. 11: Cluster G286.21+0.17 (red circle; cf. Sect. 14.4.41 and its immediate surroundings from near to far-IR. The two green circles 
mark the two bright Spitzer-Tesolved sources within it. (A close-up of the VISTA //-band image is shown in Fig. [TOl ) 




Fig. 10: Close-up of the VISTA //-band image (Fig.fTTTi around 
the two bright IRAC sources (green circles) within the cluster 
G286.21-t-0.17. 



similar case is J103754. 0-584614, which has a very faint nebu- 
lous near-lR counterpart in the VISTA image. All three are also 



detected as point-like sources in our IRAC images. They are in- 
cluded in Table IaTI (online material), but only for J103807.2- 
584511 we obtained the model parameters listed in Table IATSI 
as for the others the quality of the SED fit was not sufficient. 
Additionally, all three sources described are classified as Class I 
sources from WISE data. 

In the south-wes tern part of this extended C'^O clump, 
lYonekura etal1 ( l2005h detected a compact (r = 0.27pc) H'^CO+ 
clump (their clump No. 2); with a central density of n{H2) = 
6.8 • lO^cm"^ this is the densest of all the clumps they detected 
in their survey of the Carina Nebula complex. The Herschel- 
identified point-like source J103750. 8-584718 coincides with 
this clump. With IRAC, there are several small point-like sources 
seen to be coincident with it and an identification is therefore im- 
possible. An IRAC-detected source slightly to the south-west of 
it, J103749.3-584722, is identified as a Class I WISE cYSO. 

This area also contains the AO supergiant HD 92207. This 
star has a strong near-IR excess, was detected as a 12 \im. 
source with IRAS (IRAS 10355-5828), and is seen as a point- 
like source surrounded by nebulosity in the MIPS 24 jam image. 
Neither the star not the surrounding nebulosity can be seen in the 
Herschel images. 
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Fig. 12: Cluster G286. 38-0. 26 from near- to far-i nfrared. The large red ellipse marks the approximate outline of the C'^O clump 
No. 6 ([YAK2005] C180 6: lYonekura e'tal]|2005l) . The blue circle marks Hi^^CO+ clump No. 2 ([YAK2005] H13CO+ 2), the green 
circle indicates the shape of the arc-like nebulosity visible in the Spitzer IRAC images. 



4.4.6. Other structure 

Towards the south-west of the Gum 3 1 shell, two reflection neb- 
ulae are found, GN 10.34.5 and GN 10.31.8 (both are marked 
with white ellipses in Fig. |8]l. Both are very clearly delineated 
by cYSOs and coincide with two of the most conspicuous clus- 
ters in our field of view. GN 10.34.5 is visible in the Spitzer RGB 
image (especially the 4.5 [xm band) and seen in projection with 
a dozen cYSOs. GN 10.31.8, on the other hand, is both larger 
in angular extent and coincident with four times the number of 
cYSOs. It, too, is very conspicuous in the 4.5 i.im band. 

There also is a smaller cluster of about two dozen IRAC 
point-like sources centred at aaooo = 10:38:03, 62000 - 
-58:55:09 within the Gum 31 shell (orange circle in Fig. |8]l. 
Immediately to its w est there lies molecular cloud No. 7 of 
lYonekura et al.l (l2005h . It is accompanied by emission visible 
in the Herschel bands and a peculiarly green point-like feature 
in the IRAC RGB image, that is, strong emission in the IRAC 
4.5 i^m band. 

In the distribution of IRAC-identified cYSOs there is an- 
other notable cluster consisting of 12 candidates to the west of 
the HII region around aaooo - 10:34:27, 62000 - -58:46:45 
(orange circle in Fig. |8). This region is devoid of mid-IR emis- 
sion b ut is coincident with C'^O clump No. 1 of lYonekura et al.l 
(120051) and far-lR emission as traced by Herschel. The cYSOs 
are aligned along the western ridge of the far-lR cloud as seen 
in the Herschel image and even follow the shape of its filaments, 
broadly in the shape of an arrowhead pointing eastwards. The 
northern part is better aligned with the filament shape while the 
southern part is more randomly distributed around the filament 



itself. The border is also traced by the lYonekura et al.l (l2005l) 
'^CO intensity contours. With the WISE classification, however, 
there is nothing remarkable about that region. We find Class I 
candidate J103423.7-584531 to the north and Class 11 candidate 
J103424. 6-584749 to the south, but no appearance of clustering. 

North of the Gum 31 bubble, around a2ooo = 10:37:36, 
62000 = -58:26:36, there is a cluster of stars clearly discernible 
in the IRAC image (orange circle in Fig. |8). It lies about 2' 
to the south-west of (lYonekura et al.ll2005h C'^O clump No. 4. 
Around three dozen stars are seen within this group in pro- 
jection and it is also associated with a number of Herschel- 
identified point-like sources. Two of them, J103739. 6-582756 
and J103741. 7-582629, are part of the sample analysed in 
Sect. |5] Few stars within or around this cluster are identified as 
cYSOs with WISE. J103736.3-582655 and J103741.9-582556 
are the brightest stars in the IRAC images of the cluster and both 
identified as Class 1 candidates with WISE. There is one more 
Class 1 candidate and five Class 11 candidates distributed fairly 
evenly over the cluster 

There are several more minor clusters of ~ 5 cYSOs, no- 
tably a lways coincident with local maxima in the lYonekura et al.] 
(120051) C'^O maps. 



5. SED modelling for sources with both Herschel 
and Spitzer counterparts 

Using om Herschel and Spitzer catalogues, we were able to iden- 
tify those sources that are detected as point-like sources in both 
wavelength ranges. For Herschel we applied the restriction that 
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the sources had to be detected in at least three of the five bands, 
bringing the total number down from 9 1 sources detected in at 
least two bands to 59. This results in 16 sources overall that 
can be identified in at least three Herschel bands and at least 
one of the IRAC bands. We then compared these identifications 
to the MIPS image and performed photometry as described in 
Sect. 12.1.21 for those 6 sources where we could identify a MIPS 
counterpart. 

To extend the wavelength range of our observations, we ad- 
ditionally matched t he p oint sources analysed here with sources 
from the 2M ASS dSkrutskie et al. 2006) All-Sky Catalog of 
Point Sources (ICutri et al.H2003b . This was performed applying 
the same next-neighbour search as for the inter-band matching 
within the IRAC sources (Sect. 12.1.1b . using only sources with 
quality flags A to D. We then repeated this process with the 
WISE catalogue, where we selected only those sources that had 
a signal-to-noise ratio larger than 5. A detailed overview of all 
photometric data assembled is given in Table [ATI in the online 
material. WISE 22 |j.m and MIPS 24 |a,m fluxes show some in- 
congruity, however there is no underlying pattern as to what in 
the environs of the source may have influenced the photometry. 
This does not, however, influence the findings from the SED fits. 
In the two cases where the 22 |j.m flux appears unduly high com- 
pared to the 24 i^m flux and an SED fit is performed, leaving 
out one or the other from the fit has little or no influence on the 
best-fit model. 



5.1. Modelling of the SEDs 



For S ED-fitting we used the lonhne tooll of iRobitaille et al.l 
(I2007h . This tool compares the input observational data with 
200 000 SED models for YS Os that were precorrip uted using a 
2D radiative transfer code bv lWhitnev et al.l ( l2003l) . These mod- 
els have a wide parameter space for the properties of the central 
object and its environment. 

For the fits, the distance to all objects was assumed to be 
2.3 kpc (cf. Sect.O, and the interstellar extinction range was re- 
stricted to Ay = ... 40 mag. We assumed an uncertainty of 20% 
for all Herschel fluxes. For 2MASS, Spitzer and WISE fluxes in 
addition to the individual photometric measurement uncertainty 
as given in Table lA.ll we assumed a further systematic uncer- 
tainty of 10% due e. g. to the reliability of flux calibration. For 
IRAC, photometry varies by up to 10% due to the position of the 
point-like source within the detector array and though appropri- 
ate corrections were applied in the process, this is an additional 
source of uncertainty. 

In addition to the best-fit model, we show the range of possi- 
ble parameters that can be derived from models within the range 
of x~/N - X^sJN < 2 (with representing the number of data 
points). These models are shown as grey lines in the plots in 
Fig. [13] The resulting model parameters are listed in Table IAT21 
It gives the best-fit value together with the range constrained by 
the above criterion. The resulting SEDs are shown in Fig. [13] 
We only use fits where ;\^-/A^ for the best-fit model is smaller than 
or equal to 10.0. 

5.2. Results of SED modelling 

The results of SED fits can be highly ambiguous. Many of the 
stellar and circumstellar parameters are often poorly constrained 
because the models show a hi gh degree of degeneracy (e. g. 
iMen'shchikov & Hennindri997l) . We therefore restrict our anal- 
ysis to a few selected parameters that can be relatively well deter- 



mined from these fits. These are the total luminosity, the stellar 
mass, and the mass of the circumstellar disk and envelope. 

The best-fit masses as listed in Table lA.2] are between 1 .7 Mq 
and 6.6 Mq and even the extremes of the ranges do not exceed 
1.2 Mq to 7.1 Mq. The majority of luminosities are to be found in 
a range of ~ 100 - 300 Lq, with two notable exceptions well be- 
low that at 38 L© and 42 Lq, respectively, and one exceptionally 
luminous source at 890 Lq best-fit value, corresponding with the 
highest best-fit stellar mass in our sample. Whereas 4 of the 10 
sources sampled here exhibit best-fit envelope masses of 190 Mq 
or higher, 3 are at < 100 Mq and two lower than 50 Mq. The disk 
masses span a range of about one order of magnitude between 
~ 0.01 Mq and ~ 0.1 Mq. The highest-mass star in the sample is 
the notable exeption with a disk mass of ~ 0.001 Mq. 

In a large-scale view it is immediately noticeable that all but 
two of the sources for which Herschel counterparts to Spitzer 
point-like sources are detected are to be found within the Gum 3 1 
bubble. Although the field of our study stretches further out espe- 
cially to the west, only two sources are found outside the bubble. 
These are J103557.6-590046 and J103427.3-58461 1. 

J103806. 6-584002 is remarkable in that contrary to the vast 
majority of objects it is not located in the rim of the bubble but 
within the bubble itself, being the only specimen in our sample. 
It is placed in the very head of a pillar-like filament that extends 
from the northern rim of the bubble into it (cf. Fig. |9). With the 
methods employed in Sect. [4] we classify it as a WISE Class I 
candidate and an IRAC cYSO. In the Herschel images the fila- 
ment is rather faint, but J103806.6-584002 itself is clearly visi- 
ble as a point-like source. 



6. Sources of HH jets 



In an earlier paper (lOhlendorf et al.l2012l) we have traced a num- 
ber of Herbig-Haro jets ide ntified in the Carina Nebula com- 
plex by Smith et all ( 1201 Oal) back to their protostellar sources. 
ISmith et al.l (l2010ah also identify two HH jets and two HH jet 
candidates in the Gum 3 1 bubbl e. It should be re marked that the 
HST images on which the Sm ith et al. 1 (I2010ah study is based 
cover only a very small area of the entire Gum 31 region and 
the following analysis only represents a very small section of 
the area covered in the rest of the paper Due to the small sam- 
ple size and its very limited spatial dimensions, we cannot draw 
any conclusions about the distribution or likely number of all 
jet-driving protostars within the Gum 31 region. 

As can be seen in Fig. [14] the jet HH 1002 is almost per- 
pendicular to the ionisation front. It shows a number of dis- 
tinct features, marked bylSmith et al. (2010a) as A, B and C and 
labelled thus in our figure. Smith et al. (2010a) remark that in 
the 2MASS images they detect a very likely source located at 
0^2000 = 10:36:53.9, 62000 - -58:37:19. This matches almost ex- 
actly with the location of the object identified in IRAC images 
that is most likely to be the source of the jet: J103654.0-583720, 
which we classified as being a cYSO through its WISE and 
IRAC infrared excesses in Sect. 14.11 The Herschel observations 
trace the ionisation front well, but show no point-like source to 
be coincident with J103654. 0-583720. Therefore it could not be 
included in the SED-based analysis in Sect. [5] 

HH 1003 has a more complicated structure and ISmith et al.l 
(l20T0ah discuss it as possibly being a two-part object, made out 
of two jets in close association. Tracing back the direction of 
the bow shocks in Fig. [14] we find two point-like sources in the 
IRAC images that are very probably the sources of two different 
jets. One of them, J103652.4-583809, associated with features A 
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Fig. 13: Spectral energy distributions of those objects for which we could determine fluxes with at least three Herschel and one 
Spitzer IRAC band. Filled circles mark the input fluxes. The black line shows the best fit, and the grey lines show subsequent good 
fits. The dashed line represents the stellar photosphere corresponding to the central source of the best fitting model, as it would 
appear in the absence of circumstellar dust (but including interstellar extinction). 



and B, is identified as a cYSO from the WISE and IRAC colour- 
colour diagrams, while the other is an IRAC cYSO. Both are 
coincident with faint sources in the 2MASS images. Again, the 
Herschel image shows no point-like sources. 

The two candidate outflows, HHc-1 and HHc-2 (not to 
be confused with HHc-1 and HHc-2 within the extent of the 
central Carina Nebula which were included in the study of 
IOhlendorfet"al1l2012h . could not be traced back to any IR 
sources within the scope of our study. Following the axes of the 
jets, we cannot identify any likely emitting sources within their 
immediate surroundings. 

7. Discussion and Conclusions 

In this paper we analysed Spitzer, WISE, and Herschel data to 
investigate the cloud structure and the young stellar population 
in and around the Gum 3 1 nebula. These data provide consider- 
ably better sensitivity and spatial resolution than the previously 
available data sets. 

The Herschel far-IR maps show that the bubble surrounding 
the Gum 3 1 nebula is connected to the central parts of the Carina 
Nebula. This adds strong direct support to the assumption that 



Gum 3 1 is actually part of the Carina Nebula complex, as is also 
suggested by the matc hing C'^O radial velocities measured by 
lYonekura et all (|2005|) and other recent distance determinations 
(cf. Sect.[B. 

The bubble itself has a very sharp western edge, where 
the column density derived from Herschel measurements rises 
abruptly by at least an order of magnitude. The dust tempera- 
tures range from % 20 K in the dark clouds surrounding the H II 
region to w 30 K in the H II region and up to a; 40 K near the 
location of the O-type stars. 

The very small excess fraction seen for the mid-IR sources 
in the central stellar cluster NGC 3324 0.5%) suggests that it 
is at least several Myr old already. For the whole Gum 3 1 region, 
the Spitzer and WISE data reveal about 300/660 cYSOs. These 
objects are most likely Class I protostars or Class II sources of 
solar to intermediate mass. The 59 far-IR point-like sources re- 
vealed by the Herschel data are either pre-stellar cores or em- 
bedded (Class 0) protostars, i. e. trace a younger population of 
currently forming stars. The spatial distribution of the cYSOs is 
highly non-uniform. As we expect a contamination of around 
18% most of the widely distributed YSO population we see 
is probably due to back- and foreground stars as those would 
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Fig. 14: The Herbig-Haro jets and their probable sources as seen in four different wavelengths. From left to right: Hubble Space 
Telescope RGB image (red: WFPC2 Sll filter (673 nm), green: ACS Ha +N11 filter (658 nm), blue: WFPC2 Olll filter (502 nm); 
image credit: NASA/ESA/Hubble Heritage Team (STScl/AURA)), Spitzer IRAC RGB image with 3.6 \im in blue, 4.5 |J.m in green 
and 8.0 |J.m in red, Herschel PACS 170 yum image and VISTA //-band image. The ye llow or blue arrows indicate the broad shape and 
direction of the outflow in question, white letters mark their features as described bv lSmith et al.l(l2010al) . The probable IR-identified 
sources are marked at their Spitzer coordinates with circles in yellow or blue. 



be expected to be distributed homogeneously. Many cYSOs are 
found in rather compact clusterings, and a considerable number 
is found at the inner edge of the dusty bubble surrounding the 
H 11 region. The //erac/ieZ-identified point-like sources in partic- 
ular trace the edge of the bubble. 

This led us to assume a 'colle ct and collapse' scenari o driven 
by the O stars within the cluster (IWhitworth et al.iri994b . result- 
ing in an expected bubble size of 9 - 11 pc for an age of 1 .5 - 
2.0 Myr. This agrees very well with the ~ 10 pc diameter ob- 
served. We find four cYSOs in the very tips of small pillars ex- 
tending from the bubble rim (cf. Fig. |9]l, suggesting radiative 
triggeri ng processes very similar to what is observed in the South 
Pillars (ISmith et al.ll2010bl:lGaczkowski et al.ll2013h . 

We conclude that two different modes of triggered star- 
formation occur simultaneously in the Gum 31 region: 'collect 
and collapse', as evidenced by the bubble size and the cYSOs in 
its rim, and radiative triggering, as evidenced by cYSOs in the 
heads of pillars. 

We construct the near- to far-infrared SEDs of 17 cYSOs 
and estimate basic stellar and circumstellar parameters by com- 
parison to radiative-transfer models with good-quality fits for 10 
of them. All these cYSO s are of moderate luminosity (L % 
900 Lq), clearly suggesting that they are low- or intermediate- 
mass objects (M $ 7 Mq). This agrees with the results from our 
analysis of the cYSOs in the cent ral parts of the Carina Nebula 
complex ( iGaczkowski et al.l2013h . where we found that no high- 
mass stars are currently forming. 

We identify the driving sources of two Herbig-Haro jets 
in the western rim of the Gum 31 bubble. These sources are 
also identified as cYSOs by applying colour-selection criteria 
to IRAC or WISE photometry data. 

From the total number of cYSOs observed and the IMF 
(lKroupall2002h we can estimate a total young stellar population 
for the Gum 3 1 region. Our detection limit for cYSOs is about 



1 Mq and following the IMF there should be eight times as many 
stars below this mass (> 0. 1 M©) as above it. Correcting the 
number of cYSOs given in Sect. l4.2l for the contamination esti- 
mated there, this gives a number of ~ 5000 young stellar objects 
in the region. 

A more detailed investigation of the star formation history 
in this area requires a reliable identification of the individual 
young stars. While the infrared data presented in this paper can 
reveal protostars and young stars with circumstellar disks, most 
of the slightly older (> 2 Myr old) stars cannot be detected 
by infrared excesses. Our very recent Chandra X-ray observa- 
tions of the Gum 31 region, near-lR observations with VISTA 
and other ongoing observations will allow us to finally identify 
these stars, too. This will constitute the basis for a comprehen- 
sive multi-wavelength study of this interesting region, in a way 
similar to the recent studies of the young stel lar populations 
in the central parts of the Carina Nebula (see Townslev et al.l 
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Table A.l: Source fluxes as obtained in the Spitzer (IRAC: 3.6 i^m, 4.5 \xm, 5.8 \im, and 8.0 |^m; MIPS: 2A\xm) and Herschel (PACS; 70 and 160 fim, SPIRE: 250 jam, 
350 i_im, and 500 i_im) bands, complemented with WISE photometry from the All-Sky Data Release (3.4 i^m, 4.6 [a.m, 12 |xm, and 22 |^m) and JHKj photometry as obtained from 
the 2MASS All-Sky Catalog of Point Sources. 



\J \J LIJ. 




Fj 


Fh 




F3.4 |xm 


F3.6 \xm 


F4.5 |j.m 


F4.6 ^.m 




[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


T 1 

J LVj^Z / .0- 


-J (5401 i 








n nan _i_n ni c 




Z.jVj ±U.Ujj 


^ 1 7'3 -uO OOT 


T 1 ^^'^'^7 

J lvjjd 1 .0- 




(v^A _i_n n^Q 
O.DZ4 ±u.uoy 


1 'X -i-H 1 1 
I J.J J ±U. i i 


iZ.joU ±U.UoZ 




c ono -i-o n^7 




1 ono -i-o n-^ 1 
1 j.UuZ ±U.Uj 1 


J iUjO'4-j.Z- 


-JO J i JO 






U. /VD ±U.UZj 




1 -i-o noA 
i .JO J ±u.uzo 


1 ooA -i-O O'^n 

1 .ZZO ±U.UjU 






-jo4ZDo 


0/1 _i_n m Q 
z.Z4j ±u.uiy 


Z. 1 iZ +U.UI0 


l.jjZ ±U.UZd 






U.JJO ±U.Uj4 






coo 1 -70 
-JO J iZV 






n Q/i'j _i_n n^n 


4.0ZV ±U.U40 


1 A 7Qn -i-H n^n 

14. ly\J ±U.UjU 


A7 070 -i-O n^n 


A/1 /17 _i_n 1 Q 
04.4/ ±u. ly 


T 1 r\i.nc\r\ o 


-jojZj / 








Q -^QA -i-H nA7 


c QfiQ _i_n mn 


1 7 A1(\ -i-Ci noA 
1 /.4ZU ±U.UZO 


OA no -i-Ci 1 n 
ZO.Uz. ±U. lU 




-jo4 / jl 


u. / 1 1 +U.UZU 


Z.DlU +U.UjD 


OS 1 _i_n 


7^^ -un n77 
J. ijj ±U.U/ / 


on "^/in -i-Ci noQ 


OS T 1 -un m 1 
Zo. J lU ±u.Uj 1 


1 n ^7Q -un n7n 
lU. J /V ±u.u /u 


J IUj / ZD. /- 


-Jo4oUy 












n '^'^Q -i-O OOzL 




J103737 3- 


-584700 




0.580 ±0.021 


0.873 ±0.022 




1 925 +0 026 


6 835 +0 029 




J103739.6- 


-582756 












3.464 ±0.026 


4.411 ±0.026 


J103741.7- 


-582629 












0.659 ±0.030 




J103754.0- 


-584614 






0.876 ±0.031 


5.03 ±0.15 


9.175 ±0.029 


44.960 ±0.030 


67.63 ±0.24 


J103804.9- 


-585533 








30.11 ±0.17 


50.610 ±0.024 


158.300 ±0.033 


211.13 ±0.52 


J103806.6- 


-584002 




1.119+0.033 


4.066 ±0.039 


8.316 ±0.047 


9.277 ±0.024 


13.320 ±0.035 


12.363 ±0.059 


J103807.2- 


-584512 






8.309 ±0.065 


92.87 ±0.31 


82.700 ±0.025 


129.100 ±0.032 


186.10 ±0.46 


J103810.2- 


-584527 








7.343 ±0.073 


12.330 ±0.024 


50.940 ±0.029 


50.23 ±0.16 


J103842.1- 


-584437 




2.376 ±0.026 


12.26 ±0.074 


51.46 ±0.19 




68.850 ±0.030 


112.14 ±0.29 



Notes. The given uncertainties are the individual photometric measurement uncertainties only, as derived (for IRAC and MIPS) or as obtained from the catalogue (for 2MASS and WISE). For 
Herschel we did not obtain an uncertainty from photometry, but use an estimated total uncertainty of 20%. 



Table A. 1 : Continued. 



Source 




Fs 8 |j-m 


Fg [j.m 


F[2 yLm 


F22 |J.m 


Fi4 


F70 |j.m 


F 1 60 \Lm 


F25O ^J-ITl 




FsOO |j.m 




[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[Jy] 


[Jy] 


[Jy] 


[Jy] 


[Jy] 


J 103427. 3- 


-584611 


2.84 ±0.13 


2.21 ±0.11 


- 


188.04 ±0.42 


- 


15.8 


11.0 


10.7 


12.2 


4.47 


J103557.6- 


-590046 


- 


70.36 ±0.13 


177.29 ±0.27 


344.38 ±0.46 


- 


4.03 


7.63 


6.09 


- 


- 


J103643.2- 


-583158 


- 


- 


- 


- 


- 


- 


22.6 


14.4 


21.2 


- 


J103645.9- 


-584258 




- 


- 


~ 


- 


- 


- 


9.50 


8.54 


2.69 


J103652.4- 


-583129 


162.60 +0.12 


197.00 ±0.11 


77.32 ±0.39 


623.2 ±1.9 


755.35 ±0.14 


21.3 


13.9 


11.2 


- 


11.2 


J103700.9- 


-583237 


- 


- 


486.49 ±0.48 


1171.2 ±1.3" 


395.42 ±0.14 


- 


19.1 


19.9 


11.2 


6.96 


J103703.6- 


-584751 


35.55 ±0.15 


25.51 ±0.13 


- 


257.42 ±0.61 


181.70 ±0.13 


1.32 


6.72 


10.8 


6.32 


2.89 


J103726.7- 


-584809 


- 


- 


- 


- 


- 


- 


12.2 


11.6 


9.68 


7.10 


J103737.3- 


-584700 


17.14 ±0.14 


27.08 ±0.11 


- 


631.3 ±2.6 


149.17 ±0.14 


- 


31.4 


28.1 


29.8 


18.4 


J103739.6- 


-582756 


7.78 ±0.13 


5.47 ±0.13 


- 


164.83 ±0.38 


245.81 ±0.14 


- 


7.68 


14.1 


12.0 


5.68 


J103741.7- 


-582629 














11.1 


11.4 


9.05 


6.68 


J103754.0- 


-584614 


113.80 ±0.13 


135.80 ±0.13 


134.2 ±1.4 


1146.7 ±2.4 


858.43 ±0.14 


9.03 


35.8 


79.6 


40.7 


32.4 


J103804.9- 


-585533 


309.70 ±0.14 


364.000 ±0.099 


81.48 ±0.80 


2197.1 ±1.2 




5.91 


19.8 


11.8 


9.01 


6.23 


J103806.6- 


-584002 


17.23 ±0.12 


18.75 ±0.14 


60.29 ±0.34 


158.6 ±2.2 






10.6 


13.1 


5.70 


2.09 


J103807.2- 


-584512 


286.30 ±0.15 


526.00 ±0.14 


664.2 ±1.1 


2409.2 ±1.9 


1436.24 ±0.14 


15.2 


35.1 


26.8 






J103810.2- 


-584527 


151.40 ±0.13 


245.00 ±0.15 


276.9 ±1.1 


843.8 ±1.9 


1814.13 ±0.14 


6.32 


23.6 


43.0 


46.2 


54.2 


J103842.1- 


-584437 


111.60 ±0.11 


146.40 ±0.14 


269.07 ±0.88 


714.2 ±1.4 


505.00 ±0.13 


24.1 


14.8 


15.4 


10.2 


5.63 



Notes. The given uncertainties are the individual photometric measurement uncertainties only, as derived (for IRAC and MIPS) or as obtained from the catalogue (for 2MASS and WISE). For 
Herschel we did not obtain an uncertainty from photometry, but use an estimated total uncertainty of 20%. The extremely high 22 (xm flux compared to the 24 (j.m flux is probably due to the 
sources close proximity to a much brighter point source the contribution of which might not fully have been removed in the WISE All-Sky Data Release. 



Table A. 2; Model parameters for the sources in Gum 3 1 as obtained from the iRobitaille et al.l (l2007l) models. 



Source 


Stellar mass 








Disk mass 






Envelope mass 


Total luminosity 


Best-fit 


,,2 

X 








[Me] 










[Aiol 








[Mo] 






[Le] 




model 


N 


J 103427. 3- 


-58461 1 


5.8 


[5.8 


-7.3] 


1.02 


10" 


-1 


[4.30- 10^' 


- 1.02-10- 


] 


110 


[110 - 


170] 


320 


[320 - 


400] 


3003596 


6.9 


J103557.6- 


-590046 


1.7 


[1.2 


-2.0] 


3.91 


10- 


-2 


[1.76-10-3 


-3.91-10- 


'] 


27 


[4.2 


-60] 


42 


[27 


-68] 


3015149 


5.5 


J103645.9- 


-584258 


4.2 


[3.9 


-5.2] 


9.61 


10 


-2 


[1.09-10-3 


-9.61-10- 


'] 


46 


[18 


-46] 


160 


[140- 


190] 


3006782 


8.2 


J103726.7~ 


-584809 


6.0 


[2.6 


-6.0] 


2.64 


10" 


_2 


[1.05-10-3 


- 1.53-10- 


] 


310 


[160- 


310] 


190 


[94- 


190] 


3008699 


2.6 


J103739.6- 


-582756 


5.8 


[3.9 


-6.9] 


2.92 


10- 


2 


[2.21-10-3 


-2.16-10- 


] 


190 


[140- 


410] 


240 


[160- 


360] 


3009009 


4.8 


J103741.7- 


-582629 


3.4 


[2.6 


-6.0] 


9.53 


10- 


-2 


[1.62-10-3 


- 1.53-10- 


] 


250 


[160- 


310] 


94 


[94- 


190] 


3005296 


2.8 


J103806.6- 


-584002 


1.7 


[1.7 


-2.6] 


4.57 


10- 


-2 


[1.34-10-' 


-4.57-10-'] 


120 


[52- 


120] 


38 


[38 


-58] 


3016199 


8.8 


J103807.2- 


-584512 


6.6 


[1.7 


-7.3] 


6.77 


10- 


-4 


[6.77-10-^ 


-5.06-10- 


] 


130 


[5.3- 


310] 


890 


[220- 


890] 


3010777 


7.8 


J103842.1- 


-584437 


3.5 


[1.7 


-3.9] 


1.03 


10- 


2 


[6.45-10-3 


- 1.97-10- 


] 


220 


[13- 


220] 


170 


[130- 


190] 


3011717 


6.6 



Notes. For every model parameter the best-fit- value is given in the respective first column, followed by a range defined by the minimum and maximum value obtained from models constrained by a 
criterion. The second-to-last and last columns give the identifier of the best-fit model and hsx^lN (with A' representing the number of data points). 
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